Fidelity of WRF model in simulating heat wave events over India

The evaluation of Weather Research and Forecasting (WRF) model has been performed for simulating episodic Heat Wave (HW) events of 2015 and 2016 with varied horizontal resolutions of 27 km for the entire India (d01), 9 km for the North West (NW (d02)) and South East (SE (d03)) domain. Study compares the maximum temperature (Tmax) simulated by WRF model, using six different combination of parameterization schemes, with observations from the India Meteorological Department (IMD) during the HW events. Among the six experiments, Exp2 (i.e., combination of WSM6 microphysics (MP) together with radiation parameterization CAM, Yonsei (PBL), NOAH land surface and Grell-3D convective schemes) is found closest to the observations in reproducing the temperature. The model exhibits an uncertainty of ± 2 °C in maximum temperature (Tmax) for both the regions, suggesting regional temperature is influenced by the location and complex orography. Overall, statistical results reveal that the best performance is achieved with Exp2. Further, to understand the dynamics of rising HW intensity, two case studies of HW days along with influencing parameters like Tmax, RH and prevailing wind distribution have been simulated. Model simulated Tmax during 2015 reaches up to 44 °C in NW and SE part of India. In 2016, HW is more prevailing towards NW, while in SE region Tmax reaches upto 34–38 °C with high RH (60–85%). The comparative research made it abundantly evident that these episodic events are unique in terms of duration and geographical spread which can be used to assess the WRF performance for future projections of HW.


Description of weather research and forecasting (WRF) model
WRF is a widely used RCM, developed by the National Center for Atmospheric Research (NCAR), the National Oceanic and Atmospheric Administration (NOAA), and other communities.It is a mesoscale model, capable to simulate a wide range of atmospheric phenomena, including clouds, precipitation, winds, and temperature at a resolution ranging from few kilometers to a few hundred meters 16 .It is a numerical based model designed for both weather predictions and research, is configured into three system featuring WPS, WRF and ARW Post 17 .Present paper describes in detail how each process has been included in WRF.It incorporates both physics and dynamics component to simulate the complex interaction between the atmosphere, land surface, and ocean.WRF model have different sets of parameterization configuration for physical process hence it is suitable for a multiphysics approach [18][19][20][21] .One of the essential steps in numerical weather simulation is to choose the best sets of physics option for different weather system over various geographical region and time duration.Before long-term integration, sensitivity assessment of the model to different dynamical configurations and physical parameterization is one of the vital steps 22 .In this respect WRF model is a convenient tool as it provides plenty of schemes.In simulating the relevant climatic variables, e.g., temperature and precipitation, it is highly dependent on various physical process like; cumulus convection, planetary boundary layer (PBL), land surface, and microphysical processes [23][24][25] .Numerous studies have been undertaken with the WRF model to select the most suitable configuration for regional climate simulations [26][27][28] .WRF model capability to simulate very significant severe temperature events during HW days towards coastal regions of India.evaluated the unforeseen HW conditions over SE coastal India in May, 2015 29 .Tropospheric temperature, upper levels anticyclonic flow, higher outgoing longwave radiation, below normal precipitable water as a major synoptic feature associated with HW occurrence 30 .

General model set up and simulations
Present study utilizes Advanced Research (ARW) version 4.2 of WRF, a fully compressible and nonhydrostatic model.The WRFv4.2 offers significant advantages compared to its previous version for HW studies.It incorporates improved physical parameterizations that better represents the atmospheric process relevant to HWs, such as land surface process, boundary layer turbulence, and radiation.This leads to more accurate simulations of HW characteristics, including intensity, duration, and spatial extent.Additionally, WRFv4.2 feature enhanced computational efficiency and scalability, allowing for higher resolution simulations and the ability to capture localized HW phenomena with greater precision.The initial and lateral boundary condition for the simulation is carried out with NCEP (National Centres for Environmental Prediction) FNL reanalysis meteorological data archived at 1° × 1° grid.The boundary condition data includes the land use (USGS-24 category terrestrial; LULC-30s data), terrain elevation-2 min and vegetation fraction-10 min data.Terrain and vegetation fraction datasets are available from the standard WRF pre-processing System (WPS) v4.2.High-resolution experiments using two nested domains with horizontal resolutions of 27 km and 9 km have been used for the study is represented in Fig. 1.Three domains are selected i.e., mother domain, d01 (27 km) and daughter domains; d02 (9 km) and d03 (9 km).The detailed model descriptions and configurations are presented in Table 1.The experiment is specially performed on HW affected India regions i.e., NW (d02) and SE (d03).It provides a balanced and comprehensive approach to study the heatwave dynamics, considering their distinct climatic conditions, geographical features, and societal vulnerabilities.Domain 1 (d01) covers entire India (6-38°N and 65-98°E) and nested domain (twoway nesting) d02 extend from (24.25-31 WRF provides a number of physical configuration choices, allowing the model to be optimised for a certain weather phenomenon and location.The parameterization approaches are chosen based on earlier studies taken from different papers 31,32 .In this study six experiments considering different combination of the microphysics (MP; four schemes), planetary boundary layer (PBL; two schemes), surface layer (three schemes), shortwave and longwave radiation (RA; three schemes) and of convection (CU; two schemes) scheme are evaluated.The study implemented six distinct physical parameterization schemes, which are extensively discussed in Table 2.
In order to know its impact on the T max , RH and wind speed.initially multiple WRF simulations with different physical parameterizations have been conducted.From these sets of simulation, we have identified a reduce sets of six simulation that best represents the HWs observed in 2015.Subsequently, six simulations are conducted during HW duration of 19th-30th May 2015, with a spin-up time of 72 h.Afterward, based on robust statistical measures such as fractional bias (FB), Normalized mean square error (NMSE), Root mean square error (RMSE), correlation coefficient (r), and degree of index (d), the most optimal configuration is chosen to simulate another HW event occurring from 8th to 31st May, 2016.4][35] .Additionally, in order to examine the dynamics, various meteorological parameters have been assessed using ERA5 dataset (as shown in the supplementary section).Moreover, HW days has been identified using long term climatology of WRF simulations (2006-2014) and IMD observed (1986-2014).Study by 36 utilized a 9-year WRF climatology to predict the Indian summer monsoon during active and break periods.

Methodology
Descriptive statistics Uncertainty in model datasets have been analysed by robust statistics.It comprises the Normalized mean square error (NMSE), Root mean square error (RMSE), Fractional bias (FB), correlation coefficient (r), Index of agreement (d) and fraction of model predictions within a factor of two of observation (FA2) are explain by following equations.
where, X M : Model; X O : observations; X̅ O and X̅ M : average of respective dataset; σ X M andσ X O : standard deviation of model and observation datasets.

Taylor diagram
A Taylor diagram is a graphical tool used to assess the accuracy and reliability of climate or weather model simulations by comparing them with observation (Taylor, 2001).In essence, the Taylor diagram provides a means to assess the similarity of model simulations to observed data, taking into account various statistical metrics.By comparing the correlation coefficient and RMS differences of the model and observation, as well as their respective variances, the diagram offers insight into the accuracy and reliability of the model.
(1) where E' is the centered root-mean-square difference between the model and observation, σf and σr represents the variances of the model simulations and observation, respectively, and ρ denotes the correlation coefficient between the model and observation.Further, to identify the HW events the present study adopted following IMD criteria as shown in Table 3: During a heatwave period in India, as reported by the India Meteorological Department (IMD) or other government agencies, synoptic features typically include sustained high temperatures exceeding defined thresholds, often accompanied by the dominance of a strong anticyclone fostering subsidence and clear skies.Low humidity levels contribute to arid conditions, and stable atmospheric conditions with an inversion layer impede vertical mixing, trapping heat near the surface.Reduced wind speeds and prolonged sunshine further exacerbate the high temperatures.

Results and discussion
First, six simulations using WRF model has been conducted during HW duration of 19th-30th May, 2015 considering 72 h spin up time.The WRF model is run at 27 km with a sub-domain of 9 km.The initial and boundary conditions for the simulations is derive from the NCEP FNL dataset (http:// rda.ucar.edu/ datas ets/ ds083.2/),which has a time resolution of 6 h and a spatial resolution of approximately 1° × 1° (111 km × 111 km).

Model simulations and statistical evaluation
Spatial distribution and variability of IMD observed T max and corresponding model values during HW duration of 22nd-30th May, 2015 are shown in Fig. 2. A remarkable warming over entire India is shown with high temperature value of > 44 °C in NW, SE and central regions of around 15°N-30°N latitude.The highest temperature up to 44 °C is recorded over NW as well as SE India, whereas the lowest temperature of about 28-30 °C is recorded over the Western Himalaya and North-East region.Origination of HW and high temperature over India are likely sourced from Thar desert and Rajasthan 37 .Prevailing dry continental westerlies, lack of vegetation, subsidence and abundance of dust particles are major contributing factors for high temperature over Rajasthan and Thar Desert 38 .Further, WRF simulation for T max using different combination of parameterization schemes (as shown in Table 2) are compared with IMD observed datasets.The T max variability during HW event is well captured (24-44 °C) by Exp1, Exp2, Exp3 of and Exp5 against observed T max (28-44 °C) with slight overestimation.Model capability is however good in capturing the T max over the Indian regions except Western Himalaya.A few earlier studies, have also observed the tendency of WRF to cause lower temperatures across the Himalaya region 39 .Although the actual causes are yet unknown, they are likely linked to the cloud radiative process.However, for Exp4 and Exp6 the model is unable to perform well (− 5 °C of underestimation).This combination of configuration is not up to confidence level to consider for further study.This might be due to land surface schemes which play major role in simulating HW events.
Daily mean T max bias during 22nd-30th May 2015 is shown in Fig. 3.The bias pattern in different combination of experiments present noticeable differences over most part of the domain.It is seen that all model simulation except Exp4 and Exp6, overestimate the T max around NW and SE India by ± 1 °C and ± 2 °C, respectively.Basically, Exp2 exhibit the least biases (± 1 °C) in describing the T max against the observational data relative to the other Experiments.Thus, Exp2 which is combination of WSM6 (microphysics schemes), Yonsei (PBL schemes), NOAH (land surface scheme) with Grell-3D convective parameterization schemes outperforms other Experiments for the maximum temperature fields.The underestimation of T max (< 5 °C) in Exp4 and Exp6 over the hilly region might be due to large number of uncertainties present in climate data, which is mostly in the central and mountainous areas.Overall, Exp2 is able to simulate the spatial pattern of summer T max (24-44 °C) which is closer to the observations (28-44 °C).In general, the simulation of Exp2 shows a smaller temperature bias compared to other experiments and effectively captures the temperature range associated with HWs, as observed.Furthermore, Fig. 4g-l presents the spatial root mean square error (RMSE) statistics for the six different simulations conducted with WRF.RMSE pattern reveals the short-term performances of a model by comparison with observed datasets.A smaller value indicates the better model performance.For all the combination of Experiments, RMSE ranges between 0 °C-10 °C.Exp2 describes a best performing simulation with low RMSE values.The magnitude of RMSE over parts of NW and SE India is in order of 0-5 °C.Exp4 and Exp6 shows higher RMSE (6-10 °C) values specially in central, NW and SE regions.6i,ii it can be seen that Exp2 and Exp5 generally agree best with observations, each with almost same correlation and SD.Exp2, however, has a slightly lower RMS than Exp5 over both the stations.Further, statics of Ranchi and Rourkela stations of d03 domain have been calculated.In these stations all the configurations are clustered together, indicating similar performance for simulating the maximum temperature.Comparing all the cluster experiments, Exp4 performing slightly better than others.
Other meteorological parameters like model simulated wind speed and wind direction at 10m height (u10, v10) in d02 and d03 domain is also compared with the ground-based observation during HW period of 2015 (Fig. 7i-iv).Wind data for best combination of configuration (Exp2) is represented by wind rose plot which characterise wind speed and direction at selected stations.In d02 domain (Delhi and Jaipur) the WRF model simulate most of the prevailing strong winds in west direction, which agrees well with the ground-based observation.In d03 domain (Ranchi and Rourkela) model reproduces predominantly northwest winds.

Spatiotemporal variation of Heat wave events
A temporarily rise in maximum temperature due to deviation from normal circulation pattern leading to a HW events; it may also be accompanied by high humidity.Case studies are the most efficient method to evaluate the model performance and estimate the occurrence of extreme Heat events.Present study simulates the episodic event of HW 2015 and HW 2016 at 9 km of resolution for NW and SE Indian regions.Using IMD criteria HW days have been analysed by both IMD observed and WRF model datasets.

Episodic event of 2015 HW days
During summer 22nd to 28th May, 2015 north-western, central and eastern coastal region of India experienced severe HW condition leading loss of thousands human life in extreme high temperature condition 41 .Figure 8i,ii represents the observed and simulated temperature anomalies during the period of 11th to 31st May 2015.It clearly demonstrates the model ability to capture the high temperature over HW affected region i.e., NW and SE region.During 22nd to 30th May, HW events are identified with maximum temperature departure > 4.5 °C from the normal.Following IMD criteria, more than 4 °C temperature departure persists consecutive 3 days from IMD and WRF climatology confirm HW days over these regions.The simulated WRF and observed data  85°E.it is seen that during 2015 HW case, SE region experienced a longer (8-9 days) period of maximum temperature with deficient amount of RH over entire d03 region, which resulted in heat wave over the region.As in 42 mentioned that compared to other Indian region, south eastern India suffered a deadly 2015 heat wave with extreme temperature exceeded 45 °C in many places.As an additional diagnostic measure, the wind intensity has been calculated for both regions, revealing that the wind strength during episodes of extreme temperatures ranges from approximately 3 to 9 m/s in the d02 domain and 6 to 9 m/s in the d03 domain.Further, to understand the dynamics behind rising HW intensity during HW days (22nd-30th May, 2015) influencing parameters like T max , RH and prevailing wind distribution have been studied across entire region of India.Globally, 2015 was the warmest year on record.Model simulated multiday maximum temperature reaches up to 45 °C over NW and SE coast of India.El-Nino effects, anthropogenic climate change, late monsoon etc. are the possible reason behind those days of higher temperature.Heat events over northern and mid-latitudinal continents are strongly linked to developing El Niña and La Niña events 43 .However, low temperature ranging 26 °C-32 °C is seen towards west coast region of India.It is clearly seen that in hilly region; especially western Himalaya, model is unable to capture the data perfectly.The simulated meteorological parameters are compared with the observed ERA5 datasets and the comparison (Supplementary Fig. 3) shows good agreement (− 1 °C to + 1 °C) between the model simulated and observations.Moreover, during 2015 HW event death toll appears to be higher than in other HWs during the same period across the globe.Relative humidity is one of the influencing factors behind maximum temperature rise.Figure 10ii shows low RH over HW affected region i.e., northwest, central and south east coastal India.Over the region, low RH in order of 20-40% is seen.Due to lower relative humidity (RH) levels, the daytime temperature rises rapidly in May, as the duration of daylight is longer.Consequently, the night time temperature is unable to cool significantly.These hot days have unusual atmospheric conditions with strong, warm, gusty dry winds blowing over the NW region of India.Whereas, SE India is surrounded by oceans it experiences higher temperature generally because of hot and humid conditions.www.nature.com/scientificreports/ The combination of high outgoing longwave radiation (OLR) and the absence of convection leads to a significant increase in maximum temperature.These episodes are often intensified by a persistent dry and hot wind.The model-simulated wind distribution pattern during extreme temperature week causing HW (22nd to 30th May 2015) are presented in Fig. 10iii.Over India wind speed in range of 5 m/s is found during the HW episodes.These prevailing winds carry moisture from NW to interior regions of India.Dry air transport from the NW leading to moisture reduction along northwest and south east coast of India.It is also evident, westerly to north westerly winds majorly prevailed over most of the region of India with a cyclonic flow.Relative to SE India, NW India experienced strong and intense wind during the HW days.

Episodic event of 2016 HW
In order to further demonstrate the model reliability in predicting HW events, another case study of 2016 HW days over NW and SE regions of India at 9 km resolution are being simulated utilizing the same model configuration discussed in the preceding section.Similar to previous year (HW 2015), maximum temperature departure (12th-25th May, 2016) from IMD and WRF climatology is shown in Fig. 11i,ii,iii,iv respectively.Following IMD, HW criteria T max departure more than 4.5 °C is observed during 14th to 23rd May.It is evident from the Fig. 11 that in year 2016, NW region is more severely affected by HW event than SE Indian region.A high temperature record over India is seen on May 19, 2016.In the Phalodi city of Rajasthan, the high temperature hit 51 °C.On that day, temperatures were high in most of the Rajasthan, with most stations recording maximum temperatures exceeding 46 °C.Hovmoller plot analysing the meteorological parameters during HW 2016.The Hovmoller plot for 14th to 26th May 2016 simulated T max , RH and wind magnitude during HW days are depicted in Fig. 12.The NW (d02) region recorded the hottest day from 16th to 22nd May with a maximum temperature reaches beyond 45 °C with moderate RH of 15-60%.May 19, 2016 recorded as the hottest day across the NW India, with temperatures reaching up to 48 °C before the cooling southwest monsoon rains arrive in July.While d03 recorded the moderate temperature ranging from 35 to 44 °C with high RH reaching around 60-99%.HW in North West and central India is commonly connected by unusual blocking across the North Atlantic Ocean.In the eastern coastal region of India, heatwaves (HW) are observed to coincide with westerly anomalies over the Indian landmass.
These anomalies weaken the land-sea breeze along the coastline, as noted by Ratnam et al. 42 .Consequently, hot and dry conditions prevail due to the presence of Loo, an intense afternoon overland wind that occurs across   13i-iii.T max reaches up to 43 °C is soar towards NW regions of India, SE region is less affected in this year.There is clear humidity variation near the surface across the entire India.The RH in 2016 HW days is generally moderate towards NW region and higher towards SE region.This year appears to be normally very much humid in the summer hence the impact of most lethal HW is not solely attributed to high temperatures but also influenced by the effects of humidity.WRF simulation shows good agreement with the observed ERA5 (shown in Supplementary Fig. 4).Prevailing wind (≥ 5 m/s) with anti-clock wise wind direction is seen over major parts of the country, specially towards Bay of Bengal (BOB) and south eastern region of India.It may carry moisture from BOB region and pour heavy moisture towards eastern coastal region.

Conclusions
In this study, WRF simulations of two extreme heat wave events (i.e., May 22nd-30th, 2015 and May 11th-26th, 2016) have been conducted at d01 (27 km), d02 (9 km) and d03 (9 km) regions.The present work examines the sensitivity with different combination of configuration for WRF simulations of extreme heat events over NW, SE and entire regions of India.The initial test discarded inadequate combination of parameterizations, keeping only six valid configurations for the sensitivity analysis.The combination of WSM6 microphysics (MP) along with radiation parameterization CAM, Yonsei (PBL), NOAH land surface, and Grell-3D convective schemes proves to be effective for heat wave events over the northwest and southeast regions of India.The WSM6 microphysics scheme accurately represents the cloud and precipitation process, ensuring proper simulation of convective activity and heat wave dynamics.The radiation parameterization CAM captures the radiative transfer processes and their impact on surface heating, enabling a realistic representation of heat wave conditions.The Yonsei planetary boundary layer (PBL) scheme accurately simulates the vertical mixing of heat and moisture, crucial for capturing the atmospheric stability and temperature profiles during HWs.The NOAH land surface accounts for soil moisture and vegetation dynamics, providing accurate surface temperature and energy fluxes.Lastly, the Grell-3D convective scheme effectively simulates deep convection, which is essential for capturing the intensity and spatial distribution of convective systems associated with HW events.Together, these combinations of scheme ensure a comprehensive representation of the key processes driving HW events, making it suitable for studying HWs in the NW and SE regions of India.Based on the robust statistics best combination of simulation is selected.The T max variability during a HW event is well captured by Exp2 with a slight systematic temperature overestimation (+ 1 °C to + 2 °C), which is due to the complex interactions and uncertainties associated with soil moisture, evapotranspiration, vegetation dynamics, thermal diffusion, and model resolution.Furthermore, to understand the possible dynamics during HW days, two case studies have been analysed; namely HW 2015 (22nd-30th May) and HW 2016 (14th-23rd May), and influencing parameters such as T max , RH, and prevailing wind distribution have been studied across the entire region of India.Model-simulated multiday mean maximum temperatures during HW 2015 reach upto 44 °C in NW and SE region of India.Strong El Nino effects, anthropogenic climate change, etc. might be possible causes of the temperature rise on those days.While, the episodic occurrence of HW 2016 has also been studied, and it shows that the second to third week of May was the warmest days, with temperatures reaching up to 40 °C to 45 °C, especially in the NW part of India.Those high temperature days have unusual atmospheric conditions with strong warm, gusty dry winds blowing over the NW region of India.These two case studies considered here illustrate that deadly heat waves are complicated phenomena, in aspects of both their physical characteristics and their impacts.There are numerous meteorological reasons for HW events to occur.HW in continental interiors are often associated with persistent high-pressure systems that lead the air to stagnate.In cases like the 2015 HWs, high temperatures are experienced, the temperature and humidity parameter reveal that it is the combination of the two fields that created a rare and deadly event during 2015 and 2016 HW case study.In SE region closer to the ocean, humidity played an important role.This is especially true in the 2015 HW event as daily maximum air temperatures were high with low RH.While in HW 2016 though temperature is low but high RH may cause a HW like phenomena.As in continental interiors, extreme temperatures will continue to increase as the overall climate warms, although perhaps at a somewhat lower rate.These episodic events were clearly distinct in terms of length, intensity, and geographic dispersion, which is well reflected by WRF model according to comparative studies.Thus, using the WRF model setup to predict heat wave conditions in advance can be very helpful for future forecast and can greatly aid in managing natural disasters, effects, and disaster mitigation.

Figure 1 .
Figure 1.Selected domains for WRF simulation over India.
www.nature.com/scientificreports/Descriptivestatistics of station datasetsWith the station-based observations over India, a closer inspection of model performance is indeed mandatory to know the variability and pattern of climate behaviour.In Table4model performance for NW(d02) and SE(d03) domain have been examined using different statistics e.g., NMSE, RMSE, FB, R and d to estimate the linked biases and errors.NMSE is use to evaluate the maximum temperature and RH datasets of model.It avoids bias towards model that over predict or under predict and to emphasize the scatter in the entire datasets of sample concentrations40 .Exp2 reveals low NMSE, RMSE and FB value for temperature and RH datasets and indicates good simulation of this combination of model.For Exp2 a low NMSE value of 0.02 (Delhi), 0.03 (Jaipur) of d02 domain and 0.1(Ranchi), 0.05 (Rourkela) of d03 domain is observed.Further FB has been used to determine the overestimation or underestimation of model datasets.For the combination of Exp2, a low FB value is also seen in Delhi (0.006), Jaipur (0.0067), Ranchi (0.02) and Rourkela (0.0004) cities where model undergoes overestimation w.r.t observed datasets.For Exp2, there is a strong correlation among all the cities, with r values of 0.97, 0.95, 0.96, and 0.90 for Delhi, Jaipur, Ranchi, and Rourkela cities, respectively.Additionally, the degree of agreement between the model and observed datasets is determined by calculating the d value.Exp2 exhibits large d values (0.8-0.9), indicating a high level of agreement between both datasets.Furthermore, empirical cumulative distribution function (ECDF) plot (Fig.5) has been used to verify the simulated T max .In this plot, the data is distributed from lowest to highest values, and a best-fitted distribution line is displayed.If the fitted lines of the observed and modeled data closely align, it indicates a good fit of the data to the distribution.The ECDF values of Exp2 and Exp1 demonstrate close proximity to the observation line across all stations within the d02 (Delhi, Jaipur) and d03 (Ranchi, Rourkela) domains, as depicted in Fig.5ai-ii,biii-iv.While Exp4 and Exp6 does not match closely with observation.Further, column plot also represents the similar results such as simulated T max of Exp2 are closer to the observation over Delhi, Jaipur, Ranchi and Varanasi stations.Due to missing no. of observed data in Ranchi and Rourkela station, the data seems to be sparsely distributed (Fig.5ii).In this work simulation of WRF is crucially dependent on land surface schemes.The performance of NOAH land surface in Exp1, Exp2, Exp3 and Exp5 are overestimating the observed temperature by 2-3 °C while RUC land surface and thermal diffusion, surface clay used in Exp4 and Exp6 underestimated by − 4 to − 5 °C against the observation in NW and SE regions.The meteorological evaluation (T max ) of the WRF simulated six experiments are compared with observed datasets is represented by Taylor diagram in Fig. 6i-ii.In Fig. 6, multiple aspects of complex model are evaluated in terms of correlation, standard deviation (SD) and centred root-mean square (RMS) difference.Statistics for six experiments were computed and among all

Figure 5 .
Figure 5.The empirical cumulative distribution function (ECDF) and column plot of hourly datasets of d02 and d03 stations during heat wave events (22nd-30th May, 2015).

Figure 6 .
Figure 6.Taylor diagram represents a statistical comparison between observations and six different combinations of WRF simulated maximum temperature (°C) datasets across the Delhi, Jaipur, Ranchi, and Rourkela cities.

Figure 7 .
Figure 7.Comparison of wind roses obtained from best performing configuration of Exp2 over (i) Delhi, (ii) Jaipur, (iii) Ranchi and (iv) Rourkela with observations, during heat wave period of 22nd-30th May, 2015.

Table 1 .
Description of model configuration.

Table 2 .
Different combination of physics schemes used for maximum temperature during heat wave events.

Table 3 .
IMD criteria for Heat wave analysis.When the actual maximum temperature persists at 45 °C or above; Heat wave should be declared irrespective of the normal maximum temperature.

Table 4 .
T max (°C) and RH (%) Normalized mean square error (NMSE), root mean square error (RMSE), FB, r and d values for May, 2015 HW days for d02 and d03 selected stations.The bold values inducating the best performance of the particular experiments with low RMSE, low NMSE, low FB and high r and d values.